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Abstract A new mechanical model was developed to intro¬ 
duce the maturation process of wood cells theoretically. 
Using mechanical and physical properties of the two com¬ 
ponents of the cell wall, namely, a matrix reinforced by 
oriented cellulose microfibrils, it is possible to predict the 
relation between the anisotropic released strains and the 
microfibril angle. The model used in this study is based on 
the unified hypothesis combining the compressive stress 
generated in the cell wall matrix and the tensile stress 
originating in the cellulose microfibril as a framework. It is 
simple compared to the previously derived multilayered 
model, but it does not strictly fulfill all conditions of static 
equilibrium. Nevertheless, an excellent fit with observations 
can be obtained through varying a limited number of 
parameters. 

Key words Growth stresses • Wood cell maturation • Cell 
wall • Mechanical model ■ Fiber composites 


Introduction 

During the maturation process living cells in the differenti¬ 
ating xylem underneath the cambium are subject to biome¬ 
chanical transformations. These transformations result in a 
tendency to alter cell dimensions as the cells mature. As 
new cells are “glued” onto the older, already existing wood 
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material, these deformations are no longer free to occur, 
and maturation stresses appear in the newly formed xylem. 
As a consequence of the mechanical equilibrium, an incre¬ 
ment of internal stresses is induced in the inner part of the 
stem. Therefore, by a cumulative process, year after year, 
internal residual stresses develop inside the tree trunk. 1-4 

In this paper, we present a new mechanical model de¬ 
scribing the tendency of the wood cell to deform during the 
maturation process on the basis of the displacement 
method. The “unified hypothesis” 5-9 has been used to de¬ 
velop the model. This hypothesis takes into account both 
the “lignin swelling hypothesis” proposed by Boyd 10 and the 
“cellulose tension hypothesis” suggested by Bamber. 11 The 
mechanical models proposed in this study are then used to 
describe the relations between released growth strains and 
microfibril angles (MFAs) in the S2 layers, shown by 
Okuyama et al. and Yamamoto et al. 5-8 


Mechanical modeling of the maturation process based 
on the displacement method 

Virtual unique cell model 

According to Yamamoto’s model, a theoretical model 
of maturing wood, 9 the dimensional changes of new wood 
tissue during maturation can be more or less represented by 
the behavior of a virtual unique cell originally proposed by 
Barber. 12 In Barber’s original model, the cell is reduced to 
the S2 layer and can be considered as a long thick-walled 
circular cylinder of composite fiber-reinforced material that 
consists of the reinforcing framework of the cellulose 
microfibril (CMF) and the matrix substance of the lignin- 
hemicellulose compound. According to the interpretation 
of the reinforced-matrix hypothesis 9 (H. Yamamoto, un¬ 
published), the CMF framework exists as a bundle and the 
lignin-hemicellulose matrix exists as a skeleton in the cell 
wall; therefore, it is thought that both the CMF bundle and 
the matrix skeleton occupy the same domain in the macro¬ 
scopic limit. 
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Fig. 1. Virtual unique cell model: a long, thick-walled circular cylinder 
of composite fiber-reinforced material, a, b Cell with the microfibril 
winding in one direction (Z-helix). c Virtual unique cell model with two 
sets of oriented microfibrils winding in opposite directions around the 
cell 


The average stress tensor ctJ in the laminated material of 
the two sets of microfibril coatings (i.e., the CMF bundle) 
can be given by: 

+ a t r] ) 

Therefore, we obtain: 

°ij = C{ jkl E k i — df ( 2 ) 

As mentioned above, both the CMF bundle and the 
matrix skeleton occupy the same domain in the macro¬ 
scopic limit. Then, an approximation of the average stress 
tensor o tj in the biphasic medium (i.e., the cell wall) is given 
by 


As suggested by Barber , 12 the reinforcement in the 
model is obtained by two sets of oriented microfibrils 
winding in opposite directions around the cell (Fig. 1) with 
an angle of cp (i.e., the microfibril angle). Therefore one 
unique cell can be treated instead of a tissue of many wood 
cells. 

Constitutive equations of the cell wall material 

Mechanical presentation of the cell model 

During the maturation process each one of the two phases 
of the material - the CMF bundle and the lignin-hemicellu- 
lose matrix - shrinks or swells. Because of the perfect bond¬ 
ing between the microfibrils and the matrix, internal 
stresses are generated; and shrinkage or swelling of each 
component is no longer free to occur. If a cell was able to 
mature alone, having no interactions with the inner tissue of 
the trunk, the observable dimensional change would be 
characterized by the strain tensor e ip which is a tensorial 
unknown of the mathematical problem. 

The differences between the observable strains E :j and 
the free maturation swelling of the matrix skeleton copro¬ 
duce the mechanical strains sf = s— a” developing in the 
matrix skeleton and supposed to be elastic. The free matu¬ 
ration swelling induced in the matrix skeleton off comes 
from the mechanism of the lignin swelling hypothesis . 10 Let 
C™ kl be the elastic stiffness tensor of the matrix skeleton; 
then, in the same way as for the thermal stress phenom¬ 
enon, the associated elastic stresses become 

o’f = C&(e H - aS) (1) 

In the same way, the differences between the observable 
strains £ i; - and the free maturation swelling (or shrinkage) of 
each one of the two sets of microfibril coatings, a f(+rp) :j and 
or” jy, and the associated stresses can be related by 

/(+?) _ ("/(+?>)/ _ „/(+?)) r/Fp) = a - r /(M\ 

u ij '-ijkl \ b kt lx kl h u ij C-ijk! \b u U kl I 


Elastic properties of cell wall components 

Matrix skeleton The lignin-hemicellulose matrix as a 
skeleton is assumed to be isotropic. The components of the 
stiffness tensor can then be expressed using Lame’s param¬ 
eters X and u, which are directly connected to Young’s 
modulus E m and Poisson’s ratio v m . 

The elastic constants Qf kl of the behavior law (Eq. 1) 
are unchanged by rotation of the axes of reference; then 
nonzero terms are 

C mT = C eeee — C zzll — X + 2\i, 

f^m s-wi r-m s~im _ z>m _ o (A \ 

'-'rrde '-'rrzz QQrr Bdzz zzrr rr66 Am > V*) 

= /~' m = f~' m = .. 

^dzOz ^zrzr ^rdrd M’ 

with Lame’s parameters 

, v m E m J E m 

X = t -w-r and u = -r 

(l - 2v m )(l + v m ) 2(l + v m ) 

CMF bundle The elastic properties of the CMF bundle C\ jkh 
expressed along the main axes of the cell wall cylinder (r, 6, 
z), are more difficult to establish. The general form of elas¬ 
tic compliances (S-* u ), expressed along the main axes of each 
of the two sets of microfibril coatings, are written in Eq. (5) 
as functions of the elastic constants, assuming orthotropy. 
Then nonzero terms are 


c/* — 
0 1111 

1 /E[, 

of* 

0 2222 

1 /e{ , 53333 
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~v{ 2 /Ef 
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> ^2233 
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II 


c/* — 

°3311 
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S f * = —v f iF f 
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of* — 
°2323 

1 /gL ■ 

c f* ~ 

J 3131 

l /Gf 1 54 i 2 _ 

- 1 /g{ 2 


It has been assumed that only normal efforts along mi¬ 
crofibrils and shear efforts between them can be transmitted 
through the microfibril coating. Therefore, only E{, G s , and 
Ef il are taken as nonzero elastic parameters. Shear moduli 
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are arbitrarily estimated as functions of E{ (i.e., G ^ = G), = 
gEl where g is a constant). When we provisionally suppose 
g = 0.5, Posisson’s effects are not taken into account (i.e., 
v 12 = Wi = v 23 - v 32 = v 3 i = v i 3 = 0)- The stiffness tensor C($, 
is obtained inverting the compliance tensor S f * kl (Eq. 5). 
Under such restrictive but realistic conditions, the nonzero 
terms of (fjfa are 

C3333 = El , CI323— g'E{, C/j*3]= g ■ E{ (6) 

The components of the stiffness tensor of one of the 
microfibriFs coating Cf +q,) ijk , are first expressed along the 
main axes of the cell. They are deduced from Cf* ijk , as a 
result of the rotation of angle cp around axis 1, correspond¬ 
ing to the MFA. 

The elastic stiffness C f m , of the CMF bundle, expressed 
along the main axes of cell, is calculated from super¬ 
positioning the two sets of microfibril coatings (which elimi¬ 
nates the terms with sin cp at an odd power). 


!(+<p) _ „/F<p) = a n f ( + v) = MHO 

Orr ~ °rr Pi > a B6 a B8 

= A cos 2 cp + yS 3 sin 2 Q5 

= af ^ = jdjsin 2 <p + A cos 2 cp (11) 

<4t v) = = 2(A - A)cos 9 )sin 95 

a£+ v) = -a^ = a£ +y) = -a^ = 0 

Then, the term <5j in Eq. (2) is deduced using the formula 
<5' = (12) 

and we obtain 

d f ee = A£/anV, d L = Pi E l cos 2 9 , 

61 =6l z = d! r = dl e = 0 

The term in Eq. (2) does not depend on /3,. 


rf = 

ijkl 



+ 


ffi-v) 

r -ijkl 

2 


(7) 


An important consequence is that the axes of cylindrical 
coordinates are the main axes of the CMF bundle stiffness 
tensor (Eq. 8), whose nonzero terms are 


Crrrn E leeBl Cf zz , Cg z g z , C( rzr , Clg r g , Cfgg, 

rf C s Cf C ! Cf 

^eerrf 66zz j y ~ / zzd6'> zzrri rrzzi 


( 8 ) 


The general expressions of C{ jU are given in Appendix C, as 
functions of the elastic components C{* kh g, and cp. 


Free dimensional changes of cell wall components due to 
maturation 

According to various authors, 4 ’ 7-9 the free maturation 
swelling of the matrix as a skeleton during the maturation 
process of the cell is considered an isotropic swelling 
strain, characterized by the parameter a (>0), so the form 
of the swelling a m ij ought to be a diagonal tensor: 

diag[af j = [a, a, a} (9) 

The free maturation swelling of each microfibril coating, 
assumed to be shrinkage (<0) along the microfibril direc¬ 
tion, is characterized by a diagonal tensor, «%, when 
expressed with respect to the main axis of the microfibril 
coating. 


Classical interpretation of the mechanical problem 

Consequences of geometrical and material symmetries 

The cylindrical coordinates r, 6, z are used to express the 
various tensor components. Because of the geometry of the 
cell, the homogeneity of the material properties, and assum¬ 
ing a spatially uniform field of shrinkage (or swelling) all 
over the cell wall, this mechanical problem is independent 
of the coordinate 0, and all derivatives versus 6 are zero; 3 
I d 6 = 0. Coordinates r, 6, and z correspond to the main 
axes of material symmetry for both stiffness and maturation 
swelling tensors of the matrix (Cf kb cQ and the microfibril 
skeleton ( Cf ijkh djj). As a result of the cylindrical symmetry on 
the field of displacements, radial displacement is only a 
function of r [U r (r)|. As the twist effect is restricted by 
assumption of the axisymmetrical deformation, the tan¬ 
gential displacement is zero [U 0 (r, 6, z) = 0]. As the cell 
remains a circular cylinder after maturation, the longitudi¬ 
nal displacement is a linear function of z [U,(z) = s, ■ z], 
where s, is an unknown constant of the problem. 

In such conditions, the general expression of the strain 
components are given in Eq. (14). 

E n (r, 9, z) = E ee (r, 9, z) = — , 

v ’ dr r 

e zz (r, 6, z) = ^ = £„ (14) 

Yez{r> °> z) = Yzr{ r > A z) = Yro(r, 8, z) = 0 


diag[af) = {A, A, A} ( 10 ) 

In our particular case, the dimensional change along 
the longitudinal axis of the microfibril crystal (A) can be 
different from those in the transverse direction (A = A)- 
The free maturation swelling of respective microfibril 
coating expressed along the main axes of the cell, a§ +cr) or 
are then obtained from cf* through a rotation <p 
around axis 1. 


Introduction of Eqs. (1), (2), and (14) in Eq. (3) gives 
Eq. (15), which is the stress field in terms of displacements: 

+A—-Q+Aa (15) 

dr r 

The details of the material tensors A ir B tj , C ip and D n are 
given in Appendix B. In Eq. (15), the indices ij take only the 
values rr, 98, or zz because there are no shear stresses. 
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Formulation of the mechanical conditions 


Generally speaking, o t] must satisfy the equilibrium 
condition (Eq. 16). 


( -K> cr ee 


+ r 


do rr 
dr 


= 0 , 


da zz 

dz 


= 0 


(16) 


Introduction of Eq. (15) into Eq. (16) gives a differential 
equation of U r (Eq. 17). 


M rr + ( A _ + g rr )^L _ Bee ^ 

dr 'dr r 

A (At ~ Dee ) E « F (C m — C r fj =0 


(17) 


We can obtain a unique solution of this equation under the 
boundary condition 

M r )U. = M r )U = 0 (18) 

where R 0 and R, are the inner and outer radii of the 
cylinder, respectively. This means there is very small turgor 
pressure and no transverse restriction in the fiber. By apply¬ 
ing the solution of Eq. (17) to Eq. (15), we can determine 
the expressions of o rn o ge , and a zz , although these stresses 
contain e, as an unknown parameter. 

The maturation here is supposed to be free from the 
influence of the inner part of the trunk; therefore, no 
axial restriction is induced, and the average value of o zz is 
taken as zero. Barber 12 and Yamamoto 9 (H. Yamamoto, 
unpublished) then applied the following condition (integral 
condition) (Eq. 19) to their micromechanical analyses on 
the behavior of an isolated wood fiber. 


J o zz dA = 2njo zz rdr = 0 

cross-section Rq 


(19) 


As a consequence of the above condition, we can determine 
the expression of the unknown parameter e, (= e„). The 
components of strains e n and E ee , are deduced from Eq. (14) 
through the solution of Eq. (17) and £,. The strain of the 
transverse dimensional change of the wood fiber model e t is 
given as The expressions of s, and e, are given in 

Appendix A. 


New interpretation of the mechanical condition 

Consequence of the zero stress condition (local condition) 

Barber 12 and Yamamoto 9 (H. Yamamoto, unpublished) 
adopted the integral condition (Eq. 19) in their analysis. 
Equations (Al) and (A2) in Appendix A are based on the 
integral condition. However, in this model, the proposal is 
more restricting than the case of the integral condition. 
Directly from Eq. (15), we suppose 

°zz ( = A zz ~~ + Bzz^y - c zz + D zz e, j = 0, for v r 


This is a local condition, not an integral condition as adopted 
by Yamamoto 9 (H. Yamamoto, unpublished). The local con¬ 
dition sounds rather strange from the mechanical point of 
view. Strictly speaking, this condition is contrary to the 
equilibrium conditions of the stress field: It fulfills only a 
sufficient condition for Eq. (16) but not the necessary one. 
However, as described later, the local condition can explain 
the experimental relations between the MFA and the aniso¬ 
tropic maturation strain quite rationally. In that sense, we 
think that the local condition is a matter worthy of discussion 
for investigating the behavior of wood fiber cells during the 
maturation process. 

The general solution (Eq. 21) of the differential Eq. (20) 
gives the radial displacement 

U, = cr - B + % ~ D f l r, with U. = s,z (21) 

A zz + Hzz 

where c is an integral constant. In this case, c and e, are 
unknown parameters. We can solve these unknown 
parameters under the boundary condition (Eq. 18), and we 
obtain 


£ '(= £ zz) = 


C„{A ZZ + B zz ) 

D rr{A zz + B zl ) 


C zz (A rr + B „) 

D zz{ A rr + B rr )’ 


c = 0 


( 22 ) 


The components of strains e m e ge , are deduced from Eqs. (4) 
through (21) and (22): 




C zz - Pzz £ l 

A„ + B,. 


A 77 + B 7 


C 7 . 


D C„(A ZZ + B zz ) - C zz (A rr + B n ) 
“ D Mzz+ B zz ) - D zz {A rr + B n ) 


(23) 


These are equal to the strain of the transverse dimensional 
change of the wood fiber model, £,. It must be noted that the 
solution obtained from the local condition (Eq. 20) does not 
depend on values of the inner R 0 and outer R, cell radius. 
This is a direct consequence of the hypothesis in Eq. (20) 
(i.e., the local condition). 


Assumption on the maturation swelling of the CMF 
bundle (averaged maturation swelling condition) 

Equation (12) gives the values of d{j in Eq. (2). We propose 
the following formula to calculate them, instead of Eq. 
( 12 ): 

<5,y = C{ jkl a{, (24) 

where 


f _ a ,j a ,j 




°ij = 


(25) 


Then Eq. (2) can be expressed in the same form as for the 
matrix skeleton: 


T / = 


Cquifki a u) 


( 20 ) 
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The result is diagonal when expressed along the main axes 
of the cell: 

diagiafj) = {$, &cos 2 <p + /3 3 sin 2 cp, B 3 cos 2 w 

, ' ' (26) 
+ /3 1 sin 2 9)| 

On the basis of the new mechanical interpretation, that is, a 
combination of the local condition and the averaged matu¬ 
ration swelling condition, we calculated the behavior of the 
wood fiber model during the maturation process. 


Results and discussion 

Comparison to existing models 

The relations between cp and e h s, were calculated using the 
models proposed in this study. The results are displayed in 
Figs. 2 and 3. The values given by Yamamoto 9 were applied 
in the simulations as the parameters other than /? 1( f} 3 , and a. 
The volume ratio of the cellulose crystal and the matrix 
substance was supposed to be 40:60. The most fitting values 
of /?„ /? 3 , and a were sought to predict the experimental 
tendency by the trial and error method. The description 
given by “the unified hypothesis” with a multilayered cell 
wall is shown in Fig. 4 and compared to the experimental 
results. 

Figure 3 shows the calculated result of the model based 
on the integral condition, that is, Eq. (Al) and (A2). In this 
case, the parameters displayed in Table 1 were used, and p 


(=i? 1 //? 0 ) = 1.5 was assumed for the calculation. This model 
can explain the dependence of the transverse growth strain 
on cp when the value of g is assumed to be 0. However, the 
calculated longitudinal growth strain, which is more impor¬ 
tant for practical purposes, deviates extraordinarily from 
the experimental one except in the region of MFA smaller 
than 10 degrees. On the contrary, in the multilayered cell 
model reported by Yamamoto, 9 the calculated longitudinal 
growth strain is quantitatively consistent with the experi- 



0 20 40 60 


MFA (deg.) 

Fig. 3. Calculated released maturation strains (e h £,) and microfibril 
angle (MFA) derived from the integral condition. The most fitting 
parameters are used (Table 1). The values of g are as follows: a, 0; b, 
0.2%; c, 0.5% 



MFA (deg.) Fig. 4. Calculated released maturation strains (£,, £,) and microfibril 

angle (MFA) derived from Yamamoto’s multilayered model (unified 
Fig. 2. Calculated released maturation strains (s,, £,) and microfibril hypothesis model). Parameters used for the calculation are displayed in 
angle (MFA) derived from the local condition. The most fitting param- Table 1 (isotropic shrinkage). Shaded zones are the experimental ten- 
eters are used (Table 1). Dots are experimental results of two sugi trees dencies of two sugi trees 
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Table 1. Micromechanical parameters of the matrix and microfibrils 


Microfibril shrinkage 

E m 

(GPa) 

v" 1 

a 

(%) 

E f 

(GPa) 

A 

(%) 

ft 

(%) 

Isotropic shrinkage 

Test a 

1.2 

0.3 

0.5 

40 

a 

-0.05 

Test b 

1.2 

0.3 

0.5 

40 

a 

-0.10 

Test c 

1.2 

0.3 

0.5 

40 

a 

-0.15 

Test d 

1.2 

0.3 

0.5 

40 

a 

-0.20 

Best fit for Eqs. (Al) and (A2) 

1.2 

0.3 

0.05 

40 

a 

-0.10 

Best fit for Eqs. (22) and (23) 

1.2 

0.3 

0.05 

40 

+0.9 

-0.10 


ft has no effect on the behavior of the cell model 


mental one over a wide range of MFA, although, the 
multilayered model tends to overestimate the value of the 
transverse growth strain. 

In the case of the simulation made using Eqs. (22) and 
(23), based on the new mechanical interpretation, not only 
the longitudinal growth strain but also the transverse one 
can be predicted quantitatively as displayed in Fig. 2 by 
using the parameters shown in Table 1. In this simulation, 
we provisionally assumed g = 0.5. Also in this simulation 
the critical MFA where the longitudinal growth strain 
switches from contraction to expansion is 20-30 degrees, 
which is consistent with the experimental MFA. Thus, it has 
not been necessary to use a multilayered cell wall (e.g., the 
model reported by Yamamoto et al. 8 ' 9 ) to obtain the critical 
MFA. 

As we have seen, we can quantitatively predict the rela¬ 
tion between anisotropic growth strain and microfibril angle 
by assuming the new mechanical interpretation. Therefore, 
it is believed there is some meaning in the new mechanical 
interpretation from the viewpoint of cell wall mechanics. 

Significance of the parameter fa 

The model based on the new mechanical interpretation is 
advantageous for explaining the relations between anisotro¬ 
pic growth strains and the MFA. This advantage is derived 
not only from assuming the local condition and the aver¬ 
aged maturation swelling condition but also from introduc¬ 
ing transverse maturation swelling of the cellulose 
microfibril skeleton fa. 

Figure 5 represents the effects of fa on the generation of 
the anisotropic growth strains of the model based on the 
new mechanical interpretation. As the value of fa becomes 
larger, the curves of e, and e, tend to shift upward and 
downward, respectively, from a to d in Fig. 5. In the case of 
fa = 0 (a in Fig. 5) it is almost impossible to explain the 
longitudinal compressive stress generation in the compres¬ 
sion wood region with a large MFA. In such cases it is 
difficult to simulate the experimental tendency by searching 
fitting values of fa and a. Therefore, it is thought that intro¬ 
duction of the parameter fa is indispensable for elucidating 
the origin of the anisotropic growth stress. 

It is natural to think that fa is concerned with swelling of 
the matrix regions inside the aggregation of the cellulose 
microfibrils caused by deposition of the matrix substance. 



Fig. 5. Calculated released maturation strains (e h e,) and microfibril 
angle (MFA) derived from the local condition, in the cases supposing 
various values of ft. The most fitting parameters other than ft are 
displayed in Table 1 (best fit for Eqs. 22 and 23). Values of ft: a, 0%; ft 
0.8%; c, 1.6%; d, 2.4% 

Then, it can be presumed that not only parameter a but also 
fa are induced by maturation of the matrix substance. 10 On 
the other hand, it is thought that fa is induced in the cellu¬ 
lose microfibril crystal by the mechanism of the cellulose 
tension hypothesis. 11 From the micromechanical point of 
view, the origin of fa is independent of a and fa generation. 

Moreover, it is thought that the existence of a nonzero 
shear stiffness in each microfibril coating is indispensable 
for improving the previous models (e.g., done by 
Yamamoto et al. 7-9 Most previous models were developed 
on the basis of Barber’s model, 17 in which the microfibril 
framework has no shear stiffness. According to Barber’s 
conclusion “the shear forces across the plane are opposite in 
the two sets and canceled.” On the contrary, the existence 
of a nonzero shear stiffness in each microfibril coating (i.e., 
C[ f 23 = C{f 3] = gE{) gives nonzero shear moduli Cf znr and Cf r0re 
to the stiffness tensor of the CMF bundle in the new model, 
as shown in Appendix C. 
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Conclusion 

The new mechanical model, presented here, differs from 
previous models on three particular points: (1) The elastic 
behavior of the cellulose microfibrils bundle takes into ac¬ 
count the possibility of transmiting shear efforts between 
fibers; then the elastic shear modulus is nonzero. (2) The 
most fitting conditions governing longitudinal efforts along 
the cell are the local condition and the averaged maturation 
swelling condition. The averaged values of the stress com¬ 
ponent o zz is zero on every elementary volume, instead of 
an integral condition on the overall cross section of the cell 
in other models. (3) The new model taken into account the 
anisotropy of the elastic properties and the dimensional 
changes for both the matrix skeleton and the microfibril 
bundle. 

The model based on the new mechanical interpretation 
has three main advantages: (1) The predicted values of the 
longitudinal and transverse released strains are obtained 
by the explicit functions of the material properties of each 
component of the cell wall and the microfibril orientation. 
These functions are independent of the inner and outer 
radius of the cell cross section. (2) Under Yamamoto’s ma¬ 
terial parameters, 8,9 the critical MFA is obtained at 20-30 
degrees with a simple model. It is no longer necessary to 
introduce a multilayered model to get this result. (3) To 
obtain a perfect fit between experimental results and pre¬ 
dicted values, it is necessary to consider an anisotropic di¬ 
mensional change of the microfibrils. The new model 
suggests that the most important dimensional changes are 
located in the microfibril bundle, which shrinks along the 
fibrils and swells transversally. The swelling of the matrix 
skeleton is of less importance. The new interpretation, 
presented here, may sound strange from a mechanical view¬ 
point, but it offers several important speculations to experi¬ 
mental research on the mechanism of the unified 
hypothesis. 7 " 9 
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Appendix A: strains of the dimensional changes of the 
wood fiber model deduced by the integral condition 
(Eq. 19) 


The longitudinal strain, e,, becomes 


e *H«) = 

/ and H are composed of various variables: 


(Al) 


__ NAA^k + B„), . J ., x 
H = ■ V “ ^(p* Tl - l) 


+ 


k + 1 
N 2 {B a -AJc) ( __ 
1 - k 

r a zz + b zi 

2 (Ar — B et 


(p"* +1 - l) 


- D„) + |^ z J(p 2 - l) 


j = flK,k + B zz ) ^ +1 ^ 

P 2 (B zz ~A zz k) 


+ 


+ 


1 - k 

Az + B zz 

2 (Ar ~ 


(p~ k+1 - l) 


j(C„ -c„)~ \c zz ^p 2 - 1 ) 
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R, 


N 2 = - 
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M 
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and 


L = 


M 


De»{Ar + Bn) ' Drr{B„ + B ge ) 
At ~ Bee 

_ C rr (B ge + B n ) - C ee [A n + B„) 

Ar ~ B ee 

Transverse deformation, e t , becomes 
U r 


= \N v p k ' 1 +N 2 -p- k ^ + 


k -1 Uof, D„ 


Ar - B, 


+ | Pl .p^ + p 2 .p-^ + 


Appendix C: elastic properties of the microfibril 
skeleton 

Where c = cos<p and s = simp 

cL = c{; u = o 

C$eee = C 2222 c + ^{cLi 2CLi) c s + C 33J3 s 
= E(s 2 {s 2 + 4g-c 2 ] 

C(zzz ~ C{L C d- 2{c[ 2yi + ICyzii^C S + C 2 222 s 
= E{c 2 (c 2 + 4g-i 2 ) 

^ 2 ) CLe = Claiic — S ) + [c { 2 22 + Cl 333 — 2C223ijc S 


El 


c 2 s 2 + g-(c 2 - s 2 ) 2 


cl =ci; l3 c 2 + c{; 12 s 2 = g -E!c 2 


Appendix B: material tensors 

a = r m 4- rf r = r m + cf 

**ij K -'ijrr ' ijm u ij ij $6 ' v -' ij 661 

c„ = cm + dfj, d 9 = CL + CL 

The expressions of C m , jkl and C| H are given in Eqs. (4) and 
(8) (or Appendix C). 


CLe C( 2u c + C/i 31 5 - g-E{s 
C f n ee = C/ 9(T = C4nC + C/ 311 5 — 0 
Cqqzz = C( z $e ~ CLii^ + s ) 


+ 


(c/222 + C (333 4C 3 4 3 2 )c 2 5 2 — Elc 2 s 2 (l 4 g) 


cL = CL = clLc 2 + clLs 2 = 0 



